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Abstract 

Given a set P of n points in R'', an e-kernel K C P approximates the directional width of P 
in every direction within a relative (1 — e) factor In this paper we study the stability of e-kernels 
under dynamic insertion and deletion of points to P and by changing the approximation factor e. In 
the first case, we say an algorithm for dynamically maintaining a e-kernel is stable if at most 0(1) 
points change in K as one point is inserted or deleted from P. We describe an algorithm to main- 
tain an e-kernel of size 0{l/e'''^^^^^^) in 0(l/e^'*~^-'/^ + logn) time per update. Not only does our 
algorithm maintain a stable e-kemel, its update time is faster than any known algorithm that main- 
tains an e-kernel of size 0(l/e'^''^^^/^). Next, we show that if there is an e-kernel of P of size 
K, which may be dramatically less than 0(l/e'-'*^^-'/^), then there is an (e/2)-kernel of P of size 
0(min{l/e'-'^~^^/^, KL'i/2J log''~^(l/e)}). Moreover, there exists a point set P in R'^ and a parame- 
ter e > such that if every e-kemel of P has size at least k, then any (e/2) -kernel of P has size 
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1 Introduction 



With recent advances in sensing teciinology, massive geospatial data sets are being acquired at an unprece- 
dented rate in many application areas, including GIS, sensor networks, robotics, and spatial databases. Real- 
izing the full potential of these data sets requires developing scalable algorithms for analyzing and querying 
them. Among many interesting algorithmic developments to meet this challenge, there is an extensive 
amount of work on computing a "small summary" of large data sets that preserves certain desired properties 
of the input data and on obtaining a good trade-off between the quality of the summary and its size. A 
coreset is one example of such approximate summaries. Specifically, for an input set P and a function /, 
a coreset C C P is a subset of P (with respect to /) with the property that /(C) approximates f{P). If a 
small-size coreset C can be computed quickly (much faster than computing f{P)), then one can compute 
an approximate value of f{P) by first computing C and then computing /(C). This coreset-based approach 
has been successfully used in a wide range of geometric optimization problems over the last decade. See [2i 
for a survey. 

e-kernels. Agarwal et al. [1] introduced the notion of e-kernels and proved that it is a coreset for many 
functions. For any direction u G S'^^^, let P[u] = arg maxpgp (p, u) be the extreme point in P along 
u; uj{P, u) = {P[u] — P[—u],u) is called the directional width of P in direction u. For a given e > 0, 
K C P C R'^ is called an e-kernel of P if 

{P[u] -K[u],u) < £UJ{P,U) 

for all directions u G S"'^^]^ For simplicity, we assume e G (0, 1), because for e > 1, one can choose 
a constant number of points to form an e-kernel, and we assume d is constant. By definition, if X is an 
e-kernel of P and K is a 5-kernel of X, then K is a. {6 + e)-kernel of P. 

Agarwal et al. HI showed that there exists an e-kernel of size 0(l/e('^~^)/^) and it can be computed in 
time 0{n + l/e^*^/^). The running time was improved by Chan (E) to 0(n + l/e'^^'^/^) (see also [12J). 
In a number of applications, the input point set is being updated periodically, so algorithms have also been 
developed to maintain e-kernels dynamically. Agarwal et al. [1] had described a data structure to maintain an 
e-kernel of size ©(l/e^*^"^)/^) in (log(n)/e)*^('^) time per update. The update time was recently improved 
by Chan fj\ to 0{{l/£'-'^-'^y'^)logn + l/e'^-^/^). His approach can also maintain an e-kernel of size 
0((l/e'^) logn) with update time O(logn). If only insertions are allowed (e.g. in a streaming model), the 
size of the data structure can be improved to ©(l/e^'^-i)/^) EUa. 

In this paper we study two problems related to the stability of e-kernels: how e-kernels change as we 
update the input set or vary the value of e. 

Dynamic stability. Since the aforementioned dynamic algorithms for maintaining an e-kemel focus on 
minimizing the size of the kernel, changing a single point in the input set P may drastically change the 
resulting kernel. This is particularly undesirable when the resulting kernel is used to build a dynamic data 
structure for maintaining another information. For example, kinetic data structures (KDS) based on coresets 
have been proposed to maintain various extent measures of a set of moving points [2]. If an insertion or 
deletion of an object changes the entire summary, then one has to reconstruct the entire KDS instead of 
locally updating it. In fact, many other dynamic data structures for maintaining geometric summaries also 
suffer from this undesirable property BIIQUTTI. 

We call an e-kernel s-stable if the insertion or deletion of a point causes the e-kernel to change by at 
most s points. For brevity, if s = 0(1), we call the e-kemel to be stable. Chan's dynamic algorithm can 

^This is a slightly stronger version of the definition than defined in |T] and an e-kernel K gives a relative (1 + 2e)-approximation 
of uj{P,u) for all 11 G S''"^ (i.e. oj{K,u) < uj{P,u) < (1 + 2e)uj{K,u)). 
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be adapted to maintain a stable e-kemel of size 0((l/e"'^^) logn); see Lemma 2.6 below. An interesting 
question is whether there is an efficient algorithm for maintaining a stable e-kernel of size 0(l/e('^~^)/^), 
as points are being inserted or deleted. Maintaining a stable e-kernel dynamically is difficult for two main 
reasons. First, for an input set P, many algorithms compute e-kernels in two or more steps. They first 
construct a large e-kemel K' (e.g. see [1, 7]), and then use a more expensive algorithm to create a small 
e-kemel of K'. However, if the first algorithm is unstable, then K' may change completely each time P is 
updated. Second, all of the known e-kernel algorithms rely on first finding a "rough shape" of the input set 
P (e.g., finding a small box that contains P), estimating its fatness (5). This rough approximation is used 
crucially in the computation of the e-kemel. However, this shape is itself very unstable under insertions or 
deletions to P. Overcoming these difficulties, we prove the following in Section[2} 

Theorem 1.1. Given a parameter < e < 1, a stable e-kernel of size 0(l/e*^'^^^^/^) of a set of n points in 
can be maintained under insertions and deletions in 0(l/e('^~^^/^ + log n) time. 

Note that the update time of maintaining an e-kernel of size 

0(l/e('^"i)/2) is better than that in Q. 



Approximation stability. If the size of an e-kemel K is ©(l/e^*^"^)/^), then decreasing e changes K 
quite predictably. However, this is the worst-case bound, and it is possible that the size of K may be quite 
small, e.g., 0(1), or in general much smaller than the l/e^"^"^^/^ maximum (efficient algorithms are known 
for computing e-kernels of near-optimal size [2J). Then how much can the size increase as we reduce the 
allowable error from e to e/2? For any e > 0, let k(P, e) denote the minimum size of an e-kernel of P. 
Unlike many shape simplification problems, in which the size of simplification can change drastically as we 
reduce the value of e, we show (Section|3]l that this does not happen for e-kernels and that K{P,e/2) can be 
expressed in terms of k{P, e). 

Theorem 1.2. For any point set P and for any e > 0, 

K{P,e/2) = 0(min{K(P,e)L'^/2j iog'^-2(i/£)^ i/;.('^-i)/2})_ 

Moreover, there exist a point set P and some e > such that k(P, e/2) = 0,{k{P, e) L"^/^-! ). 



2 Dynamic Stability 

In this section we describe an algorithm that proves Theorem |1.1| The algorithm is composed of a sequence 
of modules, each with certain property. We first define the notion of anchor points and fatness of a point 
set and describe two algorithms for maintaining stable e-kernels with respect to a fixed anchor: one of them 
maintains a kemel of size 0{l/e'^~^) and the other of size 0(l/e'^'^^^^/^); the former has smaller update 
time. Next, we briefly sketch how Chan's algorithm [7| can be adapted to maintain a stable e-kemel of 
size 0{l/£^'^~^^ logn). Then we describe the algorithm for updating anchor points and maintaining a stable 
kernel as the anchors change. Finally, we put these modules together to obtain the final algorithm. We make 
the following simple observation, which will be crucial for combining different modules. 

Lemma 2.1 (Composition Lemma). IfK is an s-stable e-kernel of P and K' is an s' -stable e' -kemel of K, 
then K' is an (s • s')-stable (e + e')-kernel of P. 



Anchors and fatness of a point set. We call a point set P (3-fat if 

max uj{P,u)/uj{P,v) < (3. 

Sd — l 
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If /3 is a constant, we sometimes just say that P is fat. An arbitrary point set P can be made fat by applying 
an affine transform: we first choose a set of d + 1 anchor points A = {ao, ai, . . . , ad} using the following 
procedure of Barequet and Har-Peled (5). Choose ao arbitrarily. Let ai be the farthest point from ao- Then 
inductively, let a-i be the farthest point from the flat span(ao5 • • • , Oj-i)- (See Figure [T]) The anchor points 
A define a bounding box Ia with center at oq and orthogonal directions defined by vectors from the flat 
span(ao, • • • , flj-i) to Oj. The extents of Ia in each orthogonal direction is defined by placing each Oj on a 
bounding face and extending Ia the same distance from ao in the opposite direction. Next we perform an 
affine transform Ta on P such that the vector from the flat span(ao, . . . , ai_i) to Oj is equal to e^, where 
eo = (0, . . . , 0), ei = (1, 0, . . . , 0), . . . , ed = (0, . . . , 0, 1). This ensures that Ta{P) Q Ta{Ia) = [-1, l]'^- 
The next lemma shows that Ta{P) is fat. 

Lemma 2.2. For all u e S'^'^ and for (3d < 2'^<f/'^d\, 

uj{TA{A),u)<oj{TA{P),u)<oj{TA{lA),u)<pd-^{TA{A),u). (2.1) 

Proof. The first two inequalities follow by yl C P C I a- We can upper bound max„ (xi(Ta(/a), "u) = 
a;([— 1, 1]"^, n) < ^/d. The volume of the convex hull CQny{TA{A)) is \/d\ since it is a d-simplex and 
(T4(ao) — TA{ai), Cj) = 1 for each direction e,. We can then scale com{TA{A)) by a factor 1/2 (shrinking 
the volume by factor 1/2^^) so it fits in [0, l]'^. Now we can apply a lemma from |8| that the minimum 
width of a convex shape that is contained in [0, l]'^ is at least l/d times its d-dimensional volume, which is 
l/{2'^d\). The fatness of Ta{P) foflows from the fatness of Ta{Ia)- □ 

• : 

: • 

m m ■ • 



Figure 1. Anchor points A = {ao,ai,a2}, rectangle I a, and transform Ta applied to P; square H, two-dimensional grid 
G, and one-dimensional grid Ge on the edge e of H. 

Agarwal et al. fT] show if K is an e-kernel of P, then T{K) is an e-kernel of T{P) for any affine transform 
T, which implies that one can compute an e-kernel of T{P). We will need the following generalization of 
the definition of e-kernel. For two points sets P and Q, a subset C P is called an e-kemel of P with 
respect to Q if {P[u] - K[u],u) < eu;{Q, u) for aU u G S'^"^ 

Stable e-kernels for a fixed anchor. Let yl be a set of anchor points of P, as described above. We 
describe algorithms for maintaining stable e-kernels (with respect to yl) under the assumption that A remains 
a set of anchor points of P, i.e., A (1 P C. Ia, as P is being updated by inserting and deleting points. In 
view of the above discussion, without loss of generality, we assume I a = [—1, +1]*^ and denote it by H. 
As for the static case [1, 6], we first describe a simpler algorithm that maintains a stable e-kernel of size 
0(l/e'^~^), and then a more involved one that maintains a stable e-kemel of size 0(l/e('^~^)/^). 

Set 5 = e/^fd and draw a d-dimensional grid G inside H of size 5, i.e., the side-length of each grid cell 
is at most 6; G has 0{\/6'^) cells. For each grid cell r, let Pt = P H t. For a point j; G HI lying in a grid 
cell r, let x be the vertex of r nearest to the origin; we can view x being snapped to the vertex x. For each 
facet / of H, G induces a (d — 1) -dimensional grid G/ on /; G contains a column of cells for each cell in 
G/. For each cell A G G/, we choose (at most) one point of P as follows: let r be the nonempty grid cell 
in the column of G corresponding to A that is closest to /. We choose an arbitrary point from P,-; if there is 




3 



no nonempty cell in the column, no point is chosen. Let -^^/ be the set of chosen points. Set £ = [Jj^jj Lj. 
Agarwal et al. U proved that L is an e-kernel of P. Insertion or deletion of a point in P affects at most one 
point in Lf, and it can be updated in 0(log(l/e)) time. Hence, we obtain the following: 

Lemma 2.3. Let P be a set of n points in M'^, let A Q P be a set of anchor points of P, and Ze? < e < 1 
be a parameter. P can be preprocessed in 0{n + l/e^~^) time, so that a (2d)-stable e-kernel of P with 
respect to A of size 0{l/e'^~^) can be maintained in 0(log 1/e) time per update provided that A remains 
an anchor set of P. 

Agarwal et al. fP\ and Chan fE\ have described algorithms for computing an e-kernel of size 0(l/e('^~^)/^). 
We adapt Chan's algorithm to maintain a stable e-kernel with respect to a fixed anchor A. We begin by men- 
tioning a result of Chan that lies at the heart of his algorithm. 

Lemma 2.4 (Chan fSl). LetEen,E^<F< Efor some < t < 1, and P <^ [0 : Ef-'^ xRaset of at 
most n points. For all grid points 6 G [0 : F]'^^^ x M, the nearest neighbors of each b in P can be computed 
in time 0{n + E'^-'^F). 

We now set 7 = ^/£/c for a constant c > 1 to be used in a much sparser grid than with 5. Let C = 
[—2, +2]^^ and / be a facet of C. We draw a {d — 1) -dimensional grid on / of size 7. Assuming / lies on 
the plane Xd = —2, we choose a set Bj = {(^17, . . . , id-il, —2) G \ —[2/7] < ii, . . . , id-i < [2/7]} 
of grid points. For a subset X C P and a point b, we define iIj{X, b) = arg minx^x \\x — b\\, i.e., the point 
in X such that the snapped point is nearest to b. For a set R, iJj{X, R) = {iIj{X, r) | r G R}. There is a 
one to one mapping between the faces of C and H, so we also use / to denote the corresponding facet of H. 
Let Lf he. the set of points chosen in the previous algorithm corresponding to facet / of H for computing 
an (e/2)-kernel of P. Set Gf = ip{Lf,Bf). Chan showed that S = U/ec^/ (e/2)-kernel of L 



and thus an e-kernel of P. Scaling G and Bf appropriately and using Lemma 2.4 Gj can be computed in 
0(n + 1/e'^^'^/^) time. Hence, S can be computed in 0{n + l/e*^"^/^) time. 

Note that ip{Lf, b) can be the same for many points b G Bf, so insertion or deletion of a point in P 
(and thus in Lf) may change Gf significantly, thereby making S unstable. We circumvent this problem 
by introducing two new ideas. First, ipiLf, Bf) is computed in two stages, and second it is computed in 
an iterative manner. We describe the construction and the update algorithm for /; the same algorithm is 
repeated for all facets. 

WepartitionMinto O ( 1/7°'- 1) boxes: for J = {h. . . . ,id-i) G [-I/7, l/7]'^^inZ'^"\ we define Hj = 
[zi7, (n + 1)7] X • • • X + 1)7] X [—1, +1]. We maintain a subset X (1 Lf. Initially, we set X = 

Lf. Set Xj = X r\ Hj. We define a total order on the points of Bf. Initially, we sort Bf m lexicographic 
order, but the ordering will change as insertions and deletions are performed on P. Let (61, ... , bu) be 
the current ordering of Bf. We define a map tp : Bf ^ Lf as follows. Suppose 93(61), . . . ,(/?(6j-i) 
have been defined. Let Jj = argminj \\tl;{Xj,bi) — bi\\; here denotes the snapped point of ^(•). 
We set ip{bi) = tl:{Xj.,bi). We delete ip{bi) from X (and from Xj-) and recompute iIj{Xj-, Bf). Set 
Kf = {ip{b) \ b e b}} and % = \JfKf. Computing Ji and ip{bi) takes 0(l/e(^-^)/2) time, and, by 

Lemma : 



2.4 



7p{Xj^,Bf) can be computed in 0(|Xj| + ■ I/7) = 0{l/e^'^-^^/^) time. 

It can be proved that the map 93 and the set Kf satisfy the following properties: 

(PI) ip{bi)^ip{bj)fori^j, 

(P2) ipibi) = i;{Lf\{^ibj)\j<i},h), 

(P3) Kf^i;{Lf,Bf). 

Indeed, (PI) and (P2) follow from the construction, and (P3) follows from (P2). (P3) immediately implies 
that % is an e-kemel of P. Next, we describe the procedures for updating Kf when Lf changes. These 
procedures maintain (P1)-(P3), thereby ensuring that the algorithm maintains an e-kemel. 
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Inserting a point. Suppose a point p is inserted into Lf. We add p to X. Suppose p G HI j. We recompute 
ip{Xj, Bf). Next, we update ip{-) and % as follows. We maintain a point ^ ^ Lf. Initially, ^ is set to p. 
Suppose we have processed 61, ... , Let rj ^ Lf he. the current 92(62). If ||^ — bi\\ < \\f] — bi\\, then 
we swap ^ and ip{bi), otherwise neither ^ nor ip{bi) is updated. We then process b-i+i. After processing all 
points of Bf li ^ = p, i.e., no ip{bi) is updated, we stop. Otherwise, we add p to Kf and delete ^ from 
K/. The insertion procedure makes at most two changes in Kf, and it can be verified that (P1)-(P3) are 
maintained. 

Deleting a point. Suppose p is deleted from Lf. Suppose p G Hj. If p ^ Kf, then p ^ X. We delete p 
from X and Xj and recompute V'l-^ J) B). If p G Kf, i.e., there is a 6j G i? with p = ip{bi), then p ^ X. 
We delete p from i^j and IK, recompute ip{bi), and add the new ip{bi) to ii'j. Let G Hj; we remove 
from Xj and recompute il){Xj, Bf). We modify the ordering of Bf by moving 6i from its current 
position to the end. This is the only place where the ordering of is modified. Since 6, is now the last 
point in the ordering of Bf, the new tp{bi) does not affect any other (p{bj). The deletion procedure also 
makes at most two changes in Kf and maintains (P1)-(P3). 

Finally, insertion or deletion of a point in P causes at most one insertion plus one deletion in Lj, therefore 
we can conclude the following: 

Lemma 2.5. Let P be a set of n points in W^, A a set of anchor points of P, and < e < 1 a parameter. P 
can be preprocessed in 0{n + l/e^~^) time into a data structure so that a stable e-kernel of P with respect 
to Aofsize ©(l/e^*^"^)/^) can be maintained in 0(l/e('^~^)/^) time under insertion and deletion, provided 
that A remains an anchor set of P. 



Stabilizing Chan's dynamic algorithm. We now briefly describe how Chan's Q dynamic e-kernel 
algorithm can be adapted so that it maintains a stable e-kernel of size 0{{l/e'^~^) logn). He bypasses the 
need of fixed anchors by partitioning P into h = O(logn) layers (Pi, . . . , P^), where Pi is the inner-most 
layer and P^ is the outer-most layer, |Pj| > 'yYl'j=i+i l^il ^ constant 7 > 1, and \Ph\ = Pi is 

constructed first, and then the rest of the layers are constructed recursively with the remaining points. For 
each set Pj (for i < h) there exists a set of points Ai which serve as anchor points for Pj in the sense that they 
define a bounding box Ia^ (i.e.. Pi C lAi)- Furthermore, for all u G S'^"^ we have uj{Ai, u) < a;(Pi+i, u), 
and this remains true for a|Pi| insertions or deletions to P for a constant < a < 1. After a\Pi\ updates in 
P, the layers Pj, Pj+i, . . . , P/^ are reconstructed. Also at this point layers Pj, Pj+i, . . . , Pi-i will need to be 
reconstructed if layer Pj is scheduled to be reconstructed in fewer than a\Pi\ updates. We set Kh = Ph, and 
for i < /i, an e-kemel Ki of Pj with respect to Ai is maintained using Lemma : 
an e-kernel of P; |3<;| = 0((l/e'^-i) logn). 

When a new point p is inserted into P, it is added to the outermost layer Pj, i.e., i is the largest such value. 



2.3 



The set % = U,=i Ki is 



such that p G /a^ ■ If a point is inserted into or deleted from Pj, we update Ki using Lemma 2.3 The update 
time follows from the following lemma. 

Lemma 2.6. For any < e < 1, an e-kemel % of P of size 0((l/e'^^^) logn) can be maintained in 
0(log n) f/me, and the number of changes in % at each update is 0(1). 

Proof. Recall that if the insertion or deletion of a point does not require reconstruction of any layers, the 



update time is O(logn) and by Lemma 2.3 only 0(1) changes occur in %. We start by bounding the 
amortized time spent in reconstructing layers and their kernels. 

The kernel Ki of Pj is rebuilt if at least a|Pj| updates have taken place since the last reconstruction, if 
Ki needs to be rebuilt for some / < i, or if some Kj is rebuilt for j > i and Ki is scheduled to be rebuilt 
in fewer than a|Pj| updates to P. These conditions imply that the ith layer is rebuilt after every a\Ph\2^ 
updates where j is the smallest integer such that | Pj | < \Ph\2K The third condition acts to coordinate the 
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rebuilding of the layers so that the ith layer is not rebuilt after fewer than a\Pi\/2 updates since its last 
rebuild. 

Kh is rebuilt after a/e'^~^ updates. And \Pi\ < j^~^\Ph\. Since the entire system is rebuilt after 
0(|Pi|) = 0(n) updates, we call this interval a round. We can bound the updates to 3C in a round by 
charging 0{l/e'^~^) each time a Ki is rebuilt, which occurs at most 0(|Pi|)2/a|Pj| = @{n/\Pi\) times in 
a round. 

Thus there are 0{n) updates to the e-kernel % for every Q{n) updates to P. Thus, in an amortized sense, 
for each update to P there are 0(1) updates to %. 

This process can be de-amortized by adapting the standard techniques for de-amortizing the update time 
of a dynamic data structure [10|. If a kernel Ki is valid for k insertions or deletions to P, then we start 
construction on the next kernel Ki after k/2 insertions or deletions have taken place since the last time Ki 
was rebuilt. All insertions can be put in a queue and added to K by the time k/2 steps have transpired. All 
deletions from old Ki to new Ki are then queued and removed from K before another k/2 insertions or 
deletions. This can be done by performing 0(1) queued insertions or deletions from K each insertion or 
deletion from P. □ 



Updating anchors. We now describe the algorithm for maintaining a stable e-kernel when anchors of 
P are no longer fixed and need to be updated dynamically. Roughly speaking, we divide P into inner and 
outer subsets of points. The outer subset acts as a shield so that a stable kernel of the inner subset with 



respect to a fixed anchor can be maintained using Lemma 2.3 or 2.5 When the outer subset can no longer 



act as a shield, we reconstruct the inner and outer sets and start the algorithm again. We refer to the duration 
between two consecutive reconstruction steps as an epoch. The algorithm maintains a stable kernel within 
each epoch, and the amortized number of changes in the kernel because of reconstruction at the beginning 
of a new epoch will be 0(1). As above, we use the same de-amortization technique to make the e-kemel 
stable across epochs. We now describe the algorithm in detail. 

In the beginning of each epoch, we perform the following preprocessing. Set a = 1/10 and compute 
a CK-kemel £ of P of size O(logn) using Chan's dynamic algorithm; we do not need the stable version 
of his algorithm described above. L can be updated in O(logn) time per insertion/deletion. We choose a 
parameter m, which is set to l/e'^^^ or l/e'^^^^^/^. We create the outer subset of Pby peeling off m "layers" 
of anchor points Ai, . . . , Am- Initially, we set Pq = P. Suppose we have constructed Aq, . . . , Ai^i. Set 
P,;_i = P \ U}=i ^i' ^i^d is a-kemel of Pj-i. Next, we construct the anchor set Ai of L as described 
earlier in this section. We set Pj = Pj_i \ Ai and update L so that it is an a-kernel of Pj. Let A = [j- Ai, 
A = Am, and Pj = P \ A. Let IH = (1 + a)lA- By construction Pj C M. A forms the outer subset and 
acts as a shield for Pj, which is the inner subset. Set 6 = e/(2(l + a)(/3d)^), where /3rf is the constant in 
Lemma [Z2l 

If m = l/e'^~^ (resp. l/e^*^"^)/^), we maintain a stable (5-kemel %i of Pj with respect to A of size 0{m) 
using Lemma 2.3 (resp. Lemma 2.5 1. Set % = "Xj UA; \0C\ = 0{m). We prove below that % is an e-kernel 



of P. Let p be a point that is inserted into or deleted from P. If p G H, then we update %i using Lemma 2.3 
or 2.5 On the other hand, if p lies outside H, we insert it into or delete it from A. Once A has been updated 



m times, we end the current epoch and discard the current %. We begin a new epoch and reconstruct A, Pi, 
and %i as described above. 

The preprocessing step at the beginning of a new epoch causes 0{m) changes in % and there ai^e at least 
m updates in each epoch, therefore the algorithm maintains a stable kernel in the amortized sense. As above, 
using a de-amortization technique, we can ensure that % is stable. The correctness of the algorithm follows 
from the following lemma. 
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Lemma 2.7. % is always an e-kernel of P. 



Proof. It suffices to prove the lemma for a single epoch. Since we begin a new epoch after m updates in A, 
there is at least one i such that C C 3C. Thus we can show that Ai U %j forms an e-kemel of AiU Pj. 
For any direction u € S"^"^ 

uj{Xi,u) <Lo{Pi,u) < uj{Xj,u) + 2{l + a)6-Lo{A,u) 

< uj{%j, u) + 2{l + a)6l3a ■ uj{Pj,u) [via Lemma |2.2 1 

< uj{%i,u) + 2{l + a)5l3d-uj{Ii,u) 

< oj{%i,u) + 2{l + a)5{Pdf ■oj{Ai,u) [via Lemma [221 

< oj{%i, u) + 2(1 + a)6{^df ■ oj{Pi UAi,u). 

Thus, for any direction u G S'^^^ the extreme point of Pj U Ai is either in Pj or A.^. In the first case, %i 
approximates the width within a factor of 2(1 + a)5{f3d)'^ • uj{Pi U Ai,u) = e ■ uj{Pi U Ai, u). In the second 
case, the extreme point is in % because all of Ai is in %. Thus the set Pj U Ai has an e-kemel in %, and the 
rest of the points are also in %, so % is an e-kemel of the full set P. 

The size of % starts at 0{m) because both %i and A are of size 0{m). At most m points are inserted 
outside of I a and hence into A, thus the size of IK = %i U >l is still 0{m) after m steps. Then the epoch 
ends. □ 



Using Lemmas 2.3 and|2.5[ we can bound the update time and conclude the following. 



Lemma 2.8. For a set P of n points in M.'^ and a parameter < e < 1, there is a data structure that can 
maintain a stable e-kemel of P of size: 

(a) 0{l/£^'^~^^^'^) under insertions and deletions in time 0{n£^'^~^'^/'^ + l/e^'^~^^^'^ + logn), or 

(b) Oil/e"^-^) in time 0{n£'^-^ + logn + log(l/e)). 

Proof. We build an outer kemel of size 0{m) in 0(n + mlogn) time. It lasts for r2(m) insertions or 
deletions, so its construction time can be amortized over that many steps, and thus it costs 0{n/m + log n) 
time per insertion or deletion. 

In maintaining the inner kernel the preprocessing time can be amortized over m steps, bu t the update time 



cannot. In case (a) we maintain the inner kernel of size m = 0(1/ e^^ ^^/^) with Lemma 2.5 The update 



time is 0{n e^'^ ^^/^ + l/e*-^ ^^Z^). In case (b) we maintain the inner kemel of size m = 0(l/e'^ ^) with 



Lemma 



2.3 



The update time is 0{ne'^^^ + log(l/e)). 
The update time can be made worst case using a standard de-amortization techniques lITOl . More specif- 
ically, we start rebuilding the inner and outer kernels after m/2 steps and spread out the cost over the next 
m/2 steps. We put all of the needed insertions in a queue, inserting a constant number of points to % 
each update to P. Then after the new kernel is built, we enqueue required deletions from % and perform a 
constant number each update to P over the next m/2 steps. □ 

Putting it together. For a point set P C R'^ of size n, we can produce the best size and update 



time tradeoff for stable e-kernels by invoking Lemma 2.1 to compose three stable e-kernel algorithms. 



as illustrated in Figure |2] We first apply Lemma 2.6 to maintain a stable (e/3)- kerne l %i of P of size 



2.8 



to mamtam a sta- 



0(min{n, (1/e )logn}) with update time O(logn). We then apply Lemma 
ble (e/3)-kernel %2 of %i of size 0(l/e'^-i) with update time 0{\'Xi\e'^^^ + lo^l^Cil + log(l/e)) 



0(logn + log(l/e)). Finally we apply Lemma 2.8 again to maintain a stable (e/3)-kemel % of %2 of size 
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0(l/e(<^-i)/2) with update time 0{\%2\e^'^-^^/^ + l/e^'^-^^^ + log |3<;2|) = 0{l/e^'^-^^/'^). X is a stable 
e-kernel of P of size 0{l/e^'^~^^^'^) with update time 0(logn + l/e^'^^^^/^). This completes the proof of 
Theorem ll.il 

1^2 ^1/^'"-^)/' I X I 



p 


log n 


n 


< 



5Ci 


log n 


log(n)/e'*-i 


< 



l/e^-i ^ l/e(d-i)/2 



Figure 2. Composing stable e-kernel algorithms. 



3 Approximation Stability 



In this section we prove Theorem 1.2 We first give a short proof for the lower-bound and then a more 
involved proof of the upper bound. For the upper bound, we first develop basic ideas and prove the theorem 
in r2 and before generalizing to R'^. 



3.1 Lower Bound 

Take a cyclic polytope with n vertices and Q{n^'^/'^^ ) facets and convert it into a fat polytope T using standard 
procedures [IJ- For a parameter e > 0, we add, for each facet / of y, a point pf that is e far away from the 
facet. Let P be the set of vertices of J" together with the collection of added points. We choose e sufficiently 
small so that points in P are in convex position and all non-facet faces of "P remain as faces of conv(P). 
Then the size of an optimal e-kernel of P is at most n (by taking the vertices of "P as an e-kernel), but the 
size of an optimal e/2-kernel is at least the number of facets of P, because every point of the form has 
to be present in the kernel. The first half of the lower bound is realized with 0(1/ e^'^^^^^'^) evenly-spaced 
points on a sphere, and hence the full lower bound is proved. 



3.2 Upper Bound 

By Q, it suffices to consider the case in which P is fat and the diameter of P is normalized to 1. Let DC be 
an e-kernel of P of the smallest size. Let P = conv(D<:), and = T eB"^. We have 7 C conv(P) C 
by the definition of e-kemels. It suffices to show that there is a set IK' C P such that for 7' = conv(IK'), 
y C conv(P) C and \%'\ = 0(|3^|L°'/2J log^-s^i/^)) (iJ. 

For convenience, we assume that %' is not necessarily a subset of points in P; instead, we only require 
%' to be a subset of points in conv(P). By Caratheodory's theorem, for each point G X, we can choose a 
set Px C P of at most d + 1 points such that x G conv(Pr). We set U^eJC' desired (e/2)-kemel 

oiP' I U.e3C' Px\ < {d + l)\X'\ = 0(K(P,e)L'^/2J log^-2(l/e)). 

Initially, we add a point into %' for each point in %. If p ^ DC lies on dcoiw{P), we add p to %'. 
Otherwise we project p onto dconv{P) in a direction in which p is maximal in % and add the projected 
point to X'. Abusing the notation slightly, we use !P to denote the convex hull of these initial points. For 
simplicity, we assume CP to be a simplicial polytope. 

Decomposition of Pg \ intr T. There ai^e d types of simplices on dP. In R^ these are points and edges. 
In R"^ these are points, edges, and triangles. We can decompose \ intr P into a set of regions, each region 
(t(/) corresponding to a simplex / in CP. For each simplex / in !P let /* C S'^^^ denote the dual of / in the 
Gaussian diagram of 7. Recall that if / has dimension k (0 < k < d — 1), then /* has dimension d — l — k. 
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The region "P^ \ intr CP is partitioned into a collection of | J"! regions (where | J"! is the number of faces of all 
dimensions in 7). Each simplex / in J" corresponds to a region defined 

a{f) = {f + zu\0<z<e,uen. 

For a subsimplex t ^ f, we can similarly define a region (j(r) = {T + zu\0<z<e,u£ /*}. In R^, 
there are two types of regions: point regions and edge regions. In R^, there are three types of regions: point 
regions (see Figure [4ja)), edge regions (see Figure Qb)), and triangle regions (see Figure Qc)). 

For convenience, for any point q = q + z ■ u ^ ^if)^ where g € /, < z < e, and u ^ f* , we write 
q = q[u, z] (which intuitively reads, the point whose projection onto f is q and which is at a distance z 
above / in direction u). We also write q[v] = q + z ■ v (intuitively, q[v\ is obtained by rotating q w.r.t. / 
from direction u to direction v). Similarly, we write a simplex A[ti, z] = A © z • u, where A is a simplex 
inside /, < z < e, and u & f* , and write A[f] = K® z ■ v. 

We will proceed to prove the upper bound as follows. For each type of region cr(/) we place a bounded 
number of points from cr(/) n conv(P) into %' and then prove that all points in cr(/) n conv(P) are within 
a distance e/2 from some point in 5" = conv(3C'). We begin by introducing three ways of "gridding" 
cr{f) and then use these techniques to directly prove results for several base cases, which illustrate the main 
conceptual ideas. These base cases will already be enough to prove the results in and R'^. Finally we 
generalize this to R^ using an involved recursive construction. We set a few global values: 6 = e/12d, 
9 = 2 arcsin((^/2e), and p = 6/e. 

1: Creating layers. For a point q = q[u, z] G cr(/) we classify it depending on the value z = \q — q\. If z < 
e/2, then q is already within e/2 of We then divide the range [e/2, e] into a constant H = {e/2)/ 6 
number of cases using 11 = {hi = e/2, /i2 = e/2 + 6, . . . , hu = e — 6}. If z G [hi, /i^+i), then 
we set qh^ = q[u, hi]. We define C a{f) fl conv(P) to be the set of points that are a distance 

exactly hi from /. 

2: Discretize angles. We create a constant size 9-net Uf^h = {^^ij U2, • • •} C /* of directions with the 
following properties. (1) For each q = q[u, h] G f,h there is a direction Ui G Uf h such that the 
angle between u and is at most 9. (2) For each Uj G C/f./i there is a point =pi[ui,h] G ^'/^/i;let 
^f,h = {Pi M > 1}- Uf^h is constructed by first taking a (0/2)-net Uf of /*, then for each u'- G Uf 
choosing a point pi = qi [ui , /i] G ^' j^h where Ui is within an angle 0/2 of (if one exists), and finally 
placing Ui in Uj^h- 

3: Exponential grid. Define a set p = {do, di = {1 + p)do, . . . , dm = (1 + p)"^do} of distances where 
dm < 1 and do = 6, so m = 0(log 1/e). For a face / G y, let any r G be called a support 
point of f. Let pi, . . . ,pk be the vertices of the fc-simplex /. For each pj, and each dj G 9 (where 
di < I \pj — f\ I), let pj^i be the point at distance di from pj on the segment pjf. For each boundary 
facet F of /, define a sequence of at most m simplices Fq,Fi, . . . G conv(F U f), each a homothet of 
F, so the vertices of Fi lie on segments pjf where pj G dF (see Figure[5ja)). The translation of each 
Fi is defined so it intersects a point pj j (where pj G dF) and is as close to F as possible. This set of 
{k — 1) -simplices for each F defines the exponential grid Grj- The full grid structure is revealed as 
this is applied recursively on each Fi. 

The exponential grid Gr,A on a simplex A has two important properties for a point g G A: 

(Gl) If g G conv(F n f) lies between boundary facet F and Fq, let qo be the intersection of the line 

segment qf with Fq; then \ \q — qoW < do = 6. 
(G2) If g G conv(F n f) lies between Fi-i and Fi and the segment qf intersects Fi at qi, let qp be 

the intersection of F with the ray r^; then | \qi — q\\/\\qi — qrW < p = S/£- 
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We now describe how to handle certain simple types of regions: where / is a point or an edge. These will 
be handled the same regardless of the dimension of the problem, and they (the edge case in particular) will 
be used as important base cases for higher dimensional problems. 



Point regions. Consider a point region (j{p). For each h £ ^ create 9-net Up^h for ^p,h, so Np^^ are the 
corresponding points where each pi = p[h, m] G Np^h has m E Up^h- Put each Np^h in 

For any point q = q[u' , z] G a{p) H conv(P), let q' = q[u, h] where /i G y is the largest 

value such that h < z and u G Up^h is the closest direction to u'; set q^ = q[u', h] = q'[u']. 3 ^ 

First 1 1 — (/h 1 1 < 6 because z — h<6. Second 1 1 g/i — 1 1 < 6 because the angle between u' and \ / 

u is at most 6, and they are rotated about the point p. Thus 1 1? — f?'! | < | |g — | + | — g'| | < \ 

26 < e/2. ^ = p 

Lemma 3.1. For a point region (j{p), there exists a constant number of points Kp C <j{p) H 
conv(P) such that all points q G a{p) n conv(P) are within a distance e/2 ofconv{Kp). 

Edge regions. Consider an edge region a{e) for an edge e of T. Orient e along the x-axis. For each 
/i G ]J and u G C/e,h> let '^e,h,u be the set of points in /j within an angle 6 of u. For each ^e,h,u^ we add 
to Ke the (two) points of "^e,h,u with the largest and smallest x-coordinates, denoted by p^^ and pj^^- 




(a) a{e) in R'^ (b) top view of a{e) at height h 



Figure 3. Illustration of 2 points in X' for edge case with specific G H and u G Ue.h,e- 

For any point q = q[v, z] G a{e) n conv(P), there is a point q" = q[u, h] such that /i G is the largest 
value less than z and u £ Ue,h is the closest direction to v. Furthermore, 1 1 — 1 1 < I k — 1 1 + 1 kft — 1 1 < 
{z — h) + 2esm{9/2) = 6 + 6 = 26. We can also argue that there is a point q' = q[u', z'] G PhuPt,u' 
because if q has smaller x-coordinate than pj^^ or larger x-coordinate than p^^, then q' cannot be in '^e,h,u- 
Clearly the angle between u and u' is less than 6. This also implies that h — z' < 6. Thus H^" — (?'|| < 26, 
implying Hf? — < 45 < e/2. 

Lemma 3.2. For an edge region cr(e), there exists 0(1) points Kg C (7(e) nconv(P) such that for any point 
q = q[z, v] G a{e) Pi conv(P) there is a point p = q[h, u] G conv(i^e) such that z — h <26, \ \v — u\\ <26, 
and, in particular, \ \q — p\\ < 46 < e/2. 



3.2.1 Approximation Stability in 



For 3C C P G there are |3<C| points and edges in J". Thus combining Lemmas |3.l| and |3.2| | JC| =0(1) 
and we have proven Theorem |1.2|for d = 2. 



Theorem 3.1. For any point set P £ and for any e > Owe have k{P, e/2) / k{P, e) = 0(1). 
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3.2.2 Approximation Stability in 



Construction of %' . Now consider IK C P G R'^ and the point regions, edge regions, and triangle regions 
in the decomposition of 7^ \ intr J (see Figure|4]l. By Lemmas 3.1 and 3.2 we can add 0(|3C|) points to %' 
to account for all point and edge regions. We can now focus on the 0(|3C|) triangle regions. 




(a) / is a vertex of T (b) / is an edge of T (c) / is a facet of T 



Figure 4. Illustration of regions in the partition of 7^ \ intr T in three dimensions. 

Consider a triangle region a{t) for a triangle Hn J* (see Figure [Sj a)), t* consists of a single direction, 
the one normal to t. Let r be the highest point of a{t) n conv(i-*) in direction t* . We add r to %' and we 
create an exponential grid Gr^t with r as the support point. For each edge e E Gr^t and /i G we add the 
intersection of e[t*, h] with the boundary of a{t) n conv(P) to %' , as shown in Figurejsjb). Thus, in total 
we add 0(|3<:| log(l/e)) points to %' . 




Figure 5. Illustration to aid correctness of approximation of triangle regions in R^. 

Proof of correctness. Consider any point q = q[t* ,z\ G a{t)r\ conv(P) and associate it with a boundary 
edge e of i such that q G conv(e U f). Let q^ = q[t* , h] where /i G y is the largest height such that h < z. 
If segment qf does not intersect any edge parallel to e in Gr,t, let p = f. Otherwise, let Cj be the first 
segment parallel to e in Gr^t intersected by the ray qr, and let p be the intersection. Let p = p[t* , h] which 
must be in conv(3<;') by construction. If Cj = eo, then by (Gl) we have — p|| = ||^ — < thus 
Ik — pII < 25 < e/2 and we are done. Otherwise, let cje be the intersection of e with ray fq. By (G2) 
Hp ~ ^1 I/I |p ~ I < p = 6/e. Thus, q" = q[t* , h — ep] is below the segment qeP (see Figurejsjc)) and thus 
q" G conv(IK') since triangle ppqe is in conv(3C'). Finally, | If? — g"| | = I k — 9?t| I + I k/i — q"\ | < 25 < e/2. 



This proves Theorem 1.2 for = 3. 



Theorem 3.2. For any point set P G and for any e > Owe have k(P, e/2) / k(P, e) = 0(log 1/e). 
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3.2.3 Approximation Stability in R'^ 

Construction of . The number of regions in the decomposition of ^eVintry is |y| = 0(|3<;| L^/^J ) 
For each region cr(/), we choose a set Kf of 0(log'^^^(l/e)) points such that coiiv{Kf) C conv(P) n 
a{f) C conv(A'/) (e/2)B^. Then we set X' = U/ Kf. 



When / is a 0- or 1-simplex, we apply Lemmas |3.1| and 3.2 Otherwise, to construct Kf we create a 
recursive exponential grid on /. Specifically, for all /i G y and u G Uf^h we choose a point r = f[u, h] G 
cr(/) nconv(P) and construct an exponential grid Grj with r as the support point. Next, for each A G Grj 
we recursively construct exponential grids. That is, for all /i' G and u' G {7a, /i' we choose another support 
point r' = f'[u', h'] and construct an exponential grid G^/ a on A. At each iteration the dimension of the 
simplex in the exponential grids drops by one. We continue the recursion until we get 1-simplices. Let Qf 
be the union of all exponential grids. 

Let e be a 1-simplex in Qf. For each height /i G and direction u G Ue^h we choose two points as 
described in the construction the edge region and add them to Kf. We also place the support point of each 
A £ Qf into Kf. By construction, for a A;-simplex /, Qf contains 0(log'^^^(l/e)) simplices and thus 
\Kf \ = 0(log'=-i(l/e)). Hence \%'\ = 0{\X\W^i log'^-\l/s)). 

Proof of correctness. Let / be a A;-face of 7 (0 < k < d — 1). We need to show for any point 

q = q[u, z] G conv(P) Pi (t(/), there is a point p G conv(i^j) (specifically p'^) such that Hg — p| | < e/2. 

Before describing the technical details (mainly left to the appendix), we first provide some intuition 
regarding the proof. For any q G conv(P) n cr(/) we first consider q^ = q[u,h] where /i G is the 
largest h < z. If q^ £ conv(Kj) we are done. If qj^ ^ conv(i^j), we need to find a "helper point" q^^i 
for qk. If qk-i ^ conv(i^j) we need to recursively find a "helper point" qk-2 for qt-i, and so on until 
qm G conv(Kj). Tracing back along the recursion, we can then prove that q^ (and hence q) has a nearby 
point in conv(i^j). Note that we do not prove qj is near qj-i. Formally: 

Lemma 3.3. We can construct a sequence of helper points qj = qj[uj, hj] and simplices Aj G Qf with the 
following invariants: We can construct a sequence of helper points qj = qj[uj, hj] and simplices Aj G Qf 
with the following invariants: 

(11) qj G conv(P); 

(12) hj G and hj — /ij+i < S (for j ^ k); 

(13) qj £ Aj £ Qf and the dimension of Aj is j; and 

(14) \\uj-Uj+i\\<2e(forj^k). 

Proof. Set qk as above and A^ = /. Assume that for an index j < k, qj has been defined. Since (II) 
implies qj G ^Aj,/i 7^ there is a direction u'j_i G Uf^hj such that the angle between Uj and u'j_^ is at 
most 6. Let rj G Kf be the support point of Aj, given hj and u'j_^. Let F be the facet of Aj such that 
ijj G conv(F U fj). If the segment qjfj does not intersect any simplex in the family of (j — 1) simplices 
induced by F and rj, then let qj^i = rj and terminate the recursion (see Figure[6ja)). Otherwise, let Aj_i 
be the first such simplex intersected by the ray qjfj, and let qj^i be the intersection point (see Figure [Hj^b)). 

To determine qj^i we first find the direction Uj^i, such that qj_^ = gj_i[uj_i, z'] lies on the segment 
qjVj. Then hj-i is determined as the maximum ^ G such that h < z'. We can show (II) is satisfied 
because qj_^ must be in conv(P) because qj and rj are, and then = gj_i[nj_i, ft-j-i] lies on segment 
(ij^\(tj-\- We show (14) by \\uj — Uj-\\\ < \\uj — 'Uj_il| < Invariant (12) follows because z' > hj cosO 
and thus hj — z' < hj{l — cos 6) < e{sm6) < 6. Invariant (13) holds by construction. 

Assume the recursion does not terminate early. At j = 1, since qi G conv(P) and Ai is a line segment. 



we can apply Lemma 3.2 to cr(Ai) and find a point qo = qi[uQ, hi] G conv{Kf) such that ||gi — qo\ \ < 25 



and I |ni — no| | < 29. This completes the description of the recursive definition of helper points. □ 
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(a) 



(b) 



Figure 6. (a) Terminate early: Qj-i — Tj. (b) Recursive case: q'j-i lies on the segment qjrj. 

Let qm be the last point defined in the sequence. By construction, G conv(iirj). For each m < j < k, 
let Pj = qj [um]- We have the following key lemma, which shows that pj is close to a point e coiw{Kf). 

Lemma 3.4. For each j > m, there is a point p'- G Aj [um, hj] such that 

(1) \\p'j — PjW < j5; and 

(2) Pj = p'j - 2(j - m)5um e comr{Kf). 

Proof. We prove the lemma by induction on j. For j = m, since pm = qm ^ conv(Kj), the claim is 
trivially true by setting p'^ = pm- Assume the claim is true for some j > m. Now consider the case j + 1. 
Let Pi? be the intersection of the ray qjqj+i with df on facet F. Let Xj+i be the intersection of pFjp'j with 




Figure 7. The inductive step for proving — < {j + 1)<5 and p"+i G conv{Kf). 

Case 1: If Aj (such that qj € Aj according to (13)) is the closest facet to F, then (7^+1 lies between F 
and Aj. Thus by (Gl), we know that \ \qj+i — <ij\\ < ^- We set p'jj^i = p'j[um, ^j+i]- As such, 

\\Pj+i -Pj+i\ \ = \ \Pj - qj+i\\ < \\Pj -qjW + W^j- ^ 3^ + ^ = (i + 1)*^- 

Moreover, since /ij+i — 2{j + 1 — m)5 < hj — 2(j — m)6, Pj^i lies on the segment Pjp'j and therefore 
p'j+i G couv{Kf ). 
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Case 2: Otherwise. In this case, we have by (G2) 



7^1 



\\Pj-pF,j\\ Wqj-pF^jW ~ 

We set p'j^i = Xj+i [um, hjj^i]. First observe that 

\\p'j+i -Pj+i\ \ = Ikj+i - < \\Pj -qj\\<S. 

Furthermore, let y be the intersection of pF,jPj with pj^^xj^i (see Figure jvj^right)). We then have 

\\Pj+i - y\\ < {hj+i - hj) + \\pj -p'-W + p - \\pj -pjW <5 + 2(j - m)6 + 5 = 2(j + 1 - m)6. 
Therefore p'jj^i lies below y and as such p'jj^i G conv(Kj) since triangle p'jpFjPj G couv{Kf). □ 



We now complete the proof of Theorem 1.2 with the following lemma. 
Lemma 3.5. 1 1 <? — p';,' 1 1 < e/2. 

Proof. For q = q[u, z] and qk = q[u, h], since h < z < h + 6 then 1 1(? — Q'A:| I < f^- 

Since = then invariant (14) implies that ||nm—ii/c 1 1 < {k—m)20 < (i20 and hence < 

2esin((l/2)d26') < e2{d6/e) < 2dS. 

Finally, for j = k, Lemma 3.4 implies that | — Pk\ \ ^ k5 < d6 , that | — p'^\\ < 2{k — m)5 < 2d5, 

andp'^' e conv(i^'/). 
It follows that 

\\q-Pk\\ < \\q - QkW + \\qk - PkW + \\Pk - PkW + \\Pk - PkW 

< d + 2d6 + dS + 2d5 

< Qd6 = e/2, 

as desired. □ 



3.3 Remarks 

(1) For d = 2,3, K,{P,e/2) is only a factor of 0(1) and 0(log(l/e)), respectively, larger than e); 
therefore, the sizes of optimal e-kemels in these dimensions are relative stable. However, for d> A, 
the stability drastically reduces in the worst case because of the superlinear dependency on k{P, e). 

(2) Neither the upper nor the lower bound in the theorem is tight. For d = 3, we can prove a tighter lower 
bound of (k(P, e) log(l/(e • k{P, e)))) . We conjecture in that 

K(P,e/2) = e(K(P,e)L'^/2j \og'^-^(i/ {e'^'^-'^^/^ ■ K{P,e)))) . 
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